Quantum computing of quantum chaos and imperfection effects 
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We study numerically the imperfection effects in the quantum computing of the kicked rotator model 
in the regime of quantum chaos. It is shown that there are two types of physical characteristics: 
for one of them the quantum computation errors grow exponentially with the number of qubits in 
the computer while for the other the growth is polynomial. Certain similarity between classical and 
quantum computing errors is also discussed. 
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A great interest to quantum computers has been gen- 
erated recently by prominent theoretical results and im- 
pressive experimental progress which allowed to realise 
operations with a few qubits (see |l|] for a review). The 
most striking theoretical advantage is an enormous paral- 
lelism of quantum computing. Using the Shor algorithm 
[g| the factorization of large numbers can be done ex- 
ponentially faster on a quantum computer than by any 
known algorithm on a classical computer. Also a search 
of an item in a long list is much faster on a quantum com- 
puter as shown by Grover pi. Experimentally a variety 
of physical systems is considered for realisation of one 
qubit, viewed as a two level system, and controlled cou- 
pling between a few qubits that forms the basis for reali- 
sation of a quantum computer. These systems include ion 
traps |§,||, nuclear magnetic resonance systems Q, nu- 
clear spins with interaction controlled electronically [ 7|J|] 
or by laser pulses S, electrons in quantum dots [l0|, 
Cooper pair boxes |11|, optical lattices iQ and electrons 
floating on liquid helium fl^. As a result, a two-qubit 
gate has been experimentally realized with cold ions [pi , 
and the Grover algorithm has been performed for three 
qubits made from nuclear spins in a molecule p5|. 

It is clear that in any realistic realisation of a quan- 
tum computer a special attention should be paid to the 
imperfection effects. Indeed, the imperfections are al- 
ways present and they in principle may seriously mod- 
ify the computation results comparing to the algorithms 
based on ideal qubit operations. At present the imper- 
fection effects are tested in the numerical simulations of 
the quantum Fourier transform (QFT) [Q and the Shor 
algorithm factorization of 15 [ p6|jl7| ]. The obtained re- 
sults look to be promising for the quantum computing 
indicating that a small amount of noise does not change 
strongly the computations B even if in some cases only 
rather low level of noise is tolerable [O. However, for 
different reasons these studies do not allow to obtain an- 
alytical estimates of a tolerable imperfection level for a 
large number of qubits n. Indeed, the Shor algorithm is 
rather complicated and the capability of nowadays com- 
puters become too restrictive |19,0 . Recently the effects 



of static imperfections on the the stability of quantum 
computer hardware have been determined for a broad 
regime of parameters and it has been shown that the 
quantum hardware is sufficiently robust |18|] . However, 
these results cannot be directly generalised for a specific 
quantum algorithm operating in time. 

In the view of importance of imperfection effects we 
analyse in this paper their influence on a quantum com- 
putation of quantum chaos evolution in time. The quan- 
tum chaos in time-dependent systems was studied inten- 
sively during last two decades and it has been understood 
that the quantum interference can lead to dynamical lo- 
calisation of classical diffusive excitation in a close anal- 
ogy with the Anderson localisation in a random potential 
]19|-pi| . The study of such systems should represent a se- 
rious test for quantum computing. Indeed, in a classically 
chaotic system the numerical errors grow exponentially 
with time due to exponential local instability of motion 
which leads to chaotic diffusion in the phase space. In 
the quantum case the error growth is not so strong |2|] 
but still the dynamical localisation of quantum chaos re- 
mains rather sensitive to external perturbations and noise 
p2| , p3[ . In the view of quantum computing, our work is 
relevant to the situation when a problem is engaged with 
iterative quantum gate operations on a quantum com- 
puter, e.g. the Grover's algorithm |^. 

To investigate the imperfection effects on quantum 
computing we choose the kicked rotator model intro- 
duced in |24|. This model represents the main features 
of time-dependent quantum chaos and had been studied 
extensively in numerical simulations pftpTI and exper- 
iments with cold atoms p5| , p6[ . The unitary evolution 
operator U over the period T of the perturbation is given 

by 
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where ?i = 1 so that the commutator is [n, 0] = —i and 
the classical limit corresponds to A: — > cx), T — > while 
the classical chaos parameter K = kT remains constant. 
The operator U is given by the product of two unitary 
operators representing kick Uk = exp(— i/c cos 9) and free 



rotation JJt = exp(— iTn^/2). The dynamics is consid- 
ered on N quantum levels with periodic boundary con- 
ditions. In the classical limit the dynamics is described 
by the Chirikov standard map: 



ksme- 9^9 + Tn. 



(2) 



For K > 0.9716 the global chaos sets in with the dif- 
fusive growth n? = Dt where t is given in number of 
kicks and the diffusion rate is Z? ~ fc^/2 for K > 4.5 
pl[ . The quantum interference leads to suppression of 
this diffusion after a time scale t* Ki D and exponen- 
tial localisation of the eigenstates of U operator with the 
localisation length I « D/2 [|^,|l|. 
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The most efficient way of numerical simulation of quan- 
tum dynamics (|l|) on a classical computer is based on the 
fast Fourier transforms (FFT) between and n represen- 
tations. Indeed, the operators Uk and Ut are diagonal 
in and n representations respectively, that takes 0{N) 
multiplications for their realization. The transition be- 
tween representations is done by forward and back FFT 
with 0{N\og2N) multiplications. Thus the FFT is the 
most time consuming part in the classical computations 
of model (0). 

On the contrary the quantum computer requires only 
0(log2 A^) gate operations to perform QFT (see P 27|) 
and makes very easily the forward/back transformations 
between 9 and n representations. Hence, in this part the 
quantum computer has the exponential gain comparing 
to the classical one. However, it is not so easy to reach 
the exponential gain in the multiplication by the diagonal 
matrices Uk and Ut in 9 and n representations respec- 
tively. Of course, as for the classical computation this 



can be done in 0{N) operations. In this worst case the 
quantum computer will have 0(log2 N) gain comparing 
to the classical one. We suppose that a much better per- 
formance can be reached for the above diagonal part of 
quantum algorithm with a strong gain increase. However, 
in this paper we leave the question about maximal gain 
for future research and assume that the unitary diago- 
nal parts Uk and Ut of transformation (|I|) are performed 
by some quantum circuit exactly while imperfections are 
present only in the QFT part. Namely, for the QFT 
description in p^ (see Eqs. (14-21) there) each basic 
unitary operation Aj (one-qubit) or Bjk (two-qubit) is 
rotated on a small random angle of amplitude e <IC 1. 
At e = the operation Aj is written as no ■ (?, where 
no = (I/V2, 0, l/-\/2) and cTi's are Pauli matrices while 
with imperfection Aj — fij ■ a is achieved by choosing a 
unit vector fij tilted by an angle tj (< e) from no. For 
Bjk, we simply add a random angle of size Cjk {\^jk\ < e) 
to Ojk in Eq. (18) of |£3- These random e-rotations vary 
in time that produces an effective noise in the QFT and 
quantum computing. In this way the quantum computa- 
tion with Uq qubits models the kicked rotator dynamics 
over N = 2"« levels. 
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FIG. 1. Dependence of the second moment (n^) on time "^j- 
t for different imperfection strength e in quantum comput- ^ 

ing and different number of qubits Uq. Curves are for _p' 

Uq = 13, 12, 11 at e = 10"* from top to bottom and the 
lowest curve is for e = being the same for Uq = 13, 12, 11. 
Here k = 10, K = 5 and at £ = all probability is at n = 0. 
Inset shows the upper curve up to larger times. 
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FIG. 2. Scaling of {n^)/N'^ is shown for various values of 
e and iV = 2"": 10"* < e < 2 x 10"^ and Uq = 10 (o), 11 
(full triangle), 12 (o), 13 (full square) for fc = 10 and K = 5. 
Each point represents the averaged value over 10^ kicks for 
t < 10'* and the straight line represents the scaling given by 
(n^) — Dft ~ 5e^N^t for tq < t < t^ (see text). Here and in 
the next figures the logarithms are decimal. 

The effect of imperfections in the QFT on the second 
moment (n^), computed from the probability distribu- 
tion Wn over unperturbed levels n ((n^) = Sn"-^^n)j 
is shown in Fig. 1 for the regime of quantum chaos 
(fc > K > 1) and different number of qubits n, at 
e — 10~^ and e = 0. The data show that the noise from 
imperfections produces an effective diffusive growth of 
the second moment with the rate D^ which grows expo- 
nentially with the number of qubits. In fact the data in 



Fig. 2 show that in the regime k > K > 1 this rate is 
well described by the relation Dg « 5e^2^"'' for different 
e,nq and k ]2q |. 




FIG. 3. Probability distribution Wn over unperturbed lev- 



els for Uq — 12, k = 10 and K = 5 sd, two moments of time 
t = 100 (lower one) and t — 10^ (upper one): (a) e — 10"'* 
and (b) e = 0. Initially all probability is at n = 0. 

The physical origin of the exponential error growth in 
(n^) with Uq becomes clear from Fig. 3(a) which shows 
the probability distribution Wn at two moments of time. 
At e = the probability decays exponentially from the 
initially excited level n — due to dynamical localisation 
(see Fig. 3(b)). This decay continues up to plateau with 
level Wp ~ 10~^^ which is determined by the accuracy of 
round-off errors in the classical computer being around 
Cc ^ 10~^^. In fact these errors produce an effective 
diffusive growth so that Wp ~ ec^t and Wp is increased 
approximately by factor 10'^ when t is changed from 100 
to 10^. The classical errors have certain similarities with 
the imperfection effects in the quantum case (Fig. 3(a)). 
Indeed, the quantum errors also form a plateau Wpq at 
large n the level of which grows diffusively with time: 



Wpq oc e t. However, the quantum plateau produced 
by imperfections has certain differences comparing to 
the case of a classical computer. Namely, it is formed 
from clearly pronounced peaks located around the lev- 
els rim = ±2™ with m = 1,2, ..,71^/2. This property 



is related to the QFT structure which due to imperfec- 
tions generates transitions to Uq levels Um with proba- 
bility Wn^ (X e^. At fc ^ 1 the probability on the levels 
n 7^ "-m is much smaller than at rim (data not shown), 
but in the regime of quantum chaos with fc S> 1 each peak 
starts to take the form of the central peak with exponen- 
tial localisation as at e = 0. If the localisation length 
I is larger than the distance between nearby peaks then 
they start to overlap giving more homogeneous probabil- 
ity distribution on the plateau. In addition these peaks, 
as well as the central one at n = 0, broaden by imper- 
fection noise. Also between the initial rim peaks appear 
secondary peaks which are placed at powers of 2. 
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FIG. 4. Dependence of time scale tp on system parameters 
for 10"** < e < 2 X 10 
12 (o), 13 (full square) 

for fc = 10 and K = 5. The straight line is given by Eq. 
Inset shows scaling of normalized IPR ratio (,{(-)/£,{e — 0) 



"^, and Uq — 10 (o), 11 (full triangle), 
14 (•), 15 {+), 16 (x) and 17 (A) 



In fact QFT is performed by 0{nq^) gate operations 
7| with imperfection rotations. This imperfection noise 



oq peaks with probability Wpq 



e^riat in each 



creates riq 

peak |2|]. This leads to the diffusive growth (n^ 



N^W, 



pq 



DA with D, 



nqtl" 12. The mmrerical 
factor here is taken from the data in Fig. 2. There the 
variation of riq by 30% is too small to allow to distinguish 
Tiq prefactor in front of the exponential dependence 2^"' 
from a constant. Since at e = the second moment is 
bounded due to quantum localisation of chaos and fluc- 
tuates around (n^) « £)2 ^ 4^2 ^ ^4/4 g|i|^ the im- 



perfections strongly modify (n ) after the time scale 



D'jB^ w fc7(e''n,2^"0 



(3) 



which drops exponentially with riq. Due to the finite 
system size, the imperfection induced diffusive growth of 
(r?) is saturated around the maximal value N"^ after the 
time te K, N'^ /D^ w 2/{nq€^) (it is seen in Fig. 2 for 
large e^t). The imperfection induced diffusion exists on 
the large time interval tq <^t ^te (see Fig. 2). 



From the above estimates for Wpq growth with time it 
foUows that the probabihty in Uq peaks becomes com- 
parable with the probabihty inside the central peak 



(WpqUg ^ Uq^eH 



1) after time ij, 



l/(n,e) 



For 



t <^tp the plateau level Wpq is rather low and some char- 
acteristics should remain close to their values in absence 
of imperfections. One of them is the inverse participa- 
tion ratio (IPR) £, which is often used in the problems 
with localisation and determines the effective number of 
basis states contributing to the wave function. It is de- 
fined as ^^ Wn'^ = 1/^. By comparing the value of 
IPR in the presence of imperfections with its value at 
e = we determine the time scale tp by the condition 
£,{e)/£,{e = 0) = 1.5 for different n^, e and k. The de- 
pendence of tp on the parameters is shown in Fig. 4 for 
k = 10. According to these data 



tp « 0.33/{enq) 



(4) 



in agreement with the above estimate. Thus, the depen- 
dence of tp on Uq is polynomial. We also checked that the 
numerical coefficient C = 0.33 in (Q) does not vary signif- 
icantly with k, e.g. C « 0.32 for fc = 15 and C ~ 0.35 for 
k = 20. In our opinion this is related to the fact that the 
peaks at Wn^ are rather sparse and for large rim the dis- 
tance between them is much larger than the localization 
length / ss fc^/4. The main consequence of the relation 
(Q) is that certain physical characteristics, e.g. the IPR 
^, remain non-sensitive to imperfections in the quantum 
computing during polynomially long times. 

In conclusion, our studies of imperfection effects on 
quantum computing of the kicked rotator show that for 
certain characteristics, e.g. the second moment of the 
probability distribution, the errors grow exponentially 
with the number of qubits Uq. At the same time there are 
other characteristics, e.g. IPR ^, which are much more 
stable and for which the errors grow with Uq only poly- 
nomially. However, such stable to imperfections char- 
acteristics are essentially local and are determined only 
by a small fraction of levels of the whole Hilbert space 
iV = 2"«. In a sense the imperfections determine the pre- 
cision of quantum computations and have close similarity 
with the effect of round-off errors in the classical com- 
puter. Somehow the quantum computer with imperfec- 
tions reminds a very fast classical computer with not very 
high precision. Such a property can become rather re- 
strictive for certain computations. At the same time it is 
possible that the further development of quantum error- 
correction codes [y will allow to reach a sufficiently high 
precision of quantum computations. Also it is known 
that for the classical computer the computational efforts 
grow only polynomially with the mantissa length that 
can find certain applications for quantum computing. 
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